Effect of order-parameter fluctuations on the Halperin-Lubenski-Ma first-order 
transition in superconductors and liquid crystals 
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We show that order-parameter fluctuations in a good type-I superconductor or a liquid crystal 
always increase the size of the first-order transition. This behavior is eventually changed when 
the system crosses over to inverted- XY critical behavior, with the size of the first-order transition 
vanishing as a power law with a crossover exponent. We find good agreement between our theory 
and a recent experiment on the nematic-smectic-A first-order transition in 8CB-fOCB mixtures of 
liquid crystals. 
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More than 25 years ago, Halperin, Lubensky and Ma 
(HLM) jl) and Coleman and Weinberg [Q demonstrated 
that when a scalar field is coupled to a gauge field, the 
fluctuations of the gauge field can change the nature of 
the phase transition in the theory from continuous to 
first order. The coupling of a scalar field (or an or- 
der parameter (OP)) to a massless gauge field arises 
often in physics: the Meissner transition in supercon- 
ductors jj), the Higgs mechanism in particle physics 
and the nematic-smectic-A transition in liquid crystals 
H are among the best-known examples ||, [Q. The 
analysis of HLM initially centered on the fluctuations 
of the gauge field and neglected the OP fluctuations, 
which is justifiable for good type-I superconductors (with 
Ginzburg-Landau parameter n <C 1), where the size (to 
be precisely defined shortly) of the first-order transition 
is larger. This approximation, however, inevitably breaks 
down for strong type-II materials (k 3> 1), where neglect- 
ing the OP fluctuations would yield a first-order transi- 
tion well into the critical region. Close to four dimen- 
sions, the effect of the OP fluctuations can be studied 
using the Wilson-Fisher renormalization group 0, which, 
however, in this case leads only to "run-away" flows and 
no stable critical points. This is usually interpreted as 
a sign of a first-order transition jlj, M. Today, based 
on accumulated analytical and numerical evidence, it is 
generally believed that the transition for k 1 is again 
second order, in the so-called inverted- XY universality 
class P-|l6[. The corresponding topology of the flow of 
the coupling constants under scaling transformation is 
depicted in Fig. |l|. It may thus seem natural to expect 
that the OP fluctuations should decrease the size of the 
first-order transition, finally reducing it to zero at the 
crossover to inverted- XY critical behavior. 

The fluctuation effects in question are unfortunately 
too fine to be observable even in high-T c superconductors 
on account of the smallness of the fine-structure constant 
[f^,Q. But in liquid crystals, the coupling of the smectic 
OP to the director fluctuations is stronger, and the fluc- 
tuations at the nematic-smectic-A transition become an 
issue of central importance. A recent experiment on two- 
component liquid-crystal mixtures |19| found a surpris- 



ing result: the size of the first-order transition is larger 
than the prediction of the HLM theory. Motivated by 
this result, in this Letter we consider theoretically the 
effect of OP fluctuations on the first-order transition in 
a type-I material. We show quite generally that for a 
small enough Ginzburg-Landau parameter k, the size of 
the first-order transition is indeed always larger than the 
HLM result. Crudely, the reason is that for a good type-I 
material, the effective Ginzburg-Landau parameter k al- 
ways decreases at large scales, making the material only 
more type I and the transition more strongly first order. 
The type-I region is defined here as left of the separatrix 
in Fig. [I], where the flow is qualitatively the same as that 
of HLM near four dimensions. This behavior is eventu- 
ally changed as the separatrix between type-I (first-order) 
and type-II (inverted- XY) regimes is approached, where 
the size of the transition goes to zero as a power law, 
with a crossover exponent that characterizes the tricriti- 
cal point. Using the simple one- loop recursion relations 
for the coupling constants in an isotropic version of the 




FIG. 1. Schematic flow diagram of the quartic coupling b 
and the charge e 2 (see Eq. at the critical surface T — T* . 
Note the Gaussian (G), XY, inverted- XY (Ixy), and the tri- 
critical (T) fixed points. Good type-I materials (k < 1) lie 
far to the left of the separatrix that connects the Gaussian 
and the tricritical fixed points. 
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Ginzburg-Landau- Wilson action, we show that the ex- 
perimental data of Yethiraj and Bechhoefer are fit rather 
well by our theory. 

We are interested in a general phase transition 
described by the three-dimensional Ginzburg-Landau- 
Wilson theory for a fluctuating complex OP ^f(r) min- 
imally coupled to a fluctuating vector potential (Higgs 
scalar electrodynamics) : 



H = I d 3 f[|(V - iel(r))^(r)\ 2 + - ^ l,T " - 



|*(r)| 2 



+ ^M| 4 + ^|*(r0| 6 + \{V x A{f)f], (1) 

where a(T) = a(T — T*)/T*, and we assume the su- 
perconducting gauge, V • A = 0. For simplicity, we ne- 
glect the anisotropy inherent to liquid crystals but retain 
the sixth-order term in the Ginzburg-Landau expansion, 
which will be needed later for comparison with the ex- 
periment. Fields and lengths have been chosen so that 
the number of couplings in the theory is minimal. HLM 
showed that neglecting the fluctuations of the OP and 
integrating out the vector field A around the uniform 
OP configuration ^o, to the lowest order in charge e one 
obtains the corrected mean-field Ginzburg-Landau free- 
energy per unit volume: 



F = 



a(T) 
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in units where ksT* = 1. The negative cubic term in 
the free-energy implies a first-order transition. A use- 
ful measure of the size of the first-order transition is the 
parameter t = (T c — T*) /T* , where T c is the first-order 
transition temperature. To keep the algebra simple, we 
will assume b > and temporarily set c = in Eq. pi and 
turn c back to a finite value only later when we compare 
our results with the experiment. From Eq. ^, the size of 
the transition is then 
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where k 2 = b/2e 2 is the dimensionless Ginzburg-Landau 
parameter, and F(k) — 1/k 2 pp| . For k <C 1, the 
first-order transition thus occurs at higher temperatures, 
where the OP fluctuations may indeed be neglected. The 
HLM result becomes asymptotically correct in this limit. 

To include the OP fluctuations, we assume that the 
first-order transition in the theory (Eq. |l|) with the 
coupling constants e and b occurs at some t ^ tHLM- 
The renormalizability of the theory (Eq. implies that 
rescaling the cutoff A — > A/s by an arbitrary factor s > 1 
is equivalent to changing the coupling constants into e(s) 
and b(s) [ Bfl |, with the first-order transition for these 
renormalized couplings occurring at some new temper- 
ature to. The variation of the renormalized couplings 



e(s), b(s), and the temperature t(s) with s is described 
by differential recursion relations of the form 



d\{s) 
dln(s) 



f3 x [e 2 ( S ),b( S ),t( S )] , 



(4) 



where A = {e 2 , b, t}. If t is large enough, the OP fluctu- 
ations at the first-order transition in the rescaled theory 
indeed become negligible, and to may be approximated 
by the mean-field (HLM) expression: 
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F( K (s)). 
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Together with Eqs. [|and ^, and the boundary conditions 
t(l) = t and t(s) — to, Eq. |E| determines implicitly the 
actual size of the transition t. Clearly, the above idea is 
quite general and applicable to other weakly first-order 
transitions. First, let us demonstrate that neglecting the 
interactions between fluctuations by renormalizing the 
coupling constants only according to dimensional analy- 
sis gives just the HLM result. Power counting in Eq. |l| 
implies that t(s) — ts 2 , e 2 (s) — e 2 s, b(s) — bs. Thus, 
k(s) = k, and dividing Eqs. |^ and || gives t = tjiLM- 
More generally, for small e 2 and 6, and for t <C e 2 ,6, 
the renormalized couplings obey the differential equa- 
tions (Eq. [§ with |Q: 

(6) 



t = t(2 + ue 2 - vb + C(e 4 , b 2 , be 2 )), 



(3b = b — yb + zbe 



we 



0(e 6 ), 



0(b 3 
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,b 2 e 2 ). (8) 



The signs in the above equations are chosen so that the 
numerical coefficients u, v, x, y, z and w are positive. 
Their values in principle are non-universal and weakly 
dependent on renormalization procedure. The differen- 
tial recursion relation for the Ginzburg-Landau parame- 
ter k, for k <ti 1, may be easily obtained from Eqs. |?] and 



w 2 i /-/-)/' 4 2 2\ 



(9) 



Note that the right-hand side in the last equation is nega- 
tive, and k, if initially small enough, always decreases un- 
der renormalization. This is a consequence of the e 4 term 
in Eq. ||, which generates a negative quartic coupling and 
tends to drive the transition first order. The same term 
is responsible for the runaway flows near four dimensions 
in the original HLM analysis. For small charge e, one 
may neglect the nonlinear terms and solve Eqs. |^ and 
to find, for small k, 



t 



tHLM 



(10) 
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Since k(s) < k, the solution to the above equation al- 
ways gives t > tjiLM- I n sum, while near the Gaussian 
fixed point t and e 2 approximately scale according to 
dimensional analysis, a small Ginzburg-Landau parame- 
ter k acquires a small negative dimension and renormal- 
izes downwards, thus increasing the size of the first-order 
transition. 

Assuming the general topology of the flow diagram as 
shown in Fig. p] the size of the first-order transition will 
go to zero as the coupling constants approach the separa- 
trix between type-I and type-11 regimes. If the separatrix 
at the critical surface t = lies at some b c (e), close to 
the separatrix one finds 

toe [6 c (e)-6] 1/0 , (11) 

where 4> = rv is the crossover exponent, with v^ 1 = 
dPt/dt, taken at the tricritical fixed point. Here, r is the 
(positive) scaling dimension of the second relevant scal- 
ing variable at the tricritical point, which in general is a 
linear combination of e 2 , 6, and t. Near four dimensions, 
our result agrees with that of HLM: the transition is then 
always first order, unless e = 0. The role of the tricriti- 
cal point near four dimensions is played by the XY fixed 
point, which has two relevant directions, t and e 2 , with 
v = v xy and r = e, where e = 4 — d. Instead of Eq. [Il], 
one then has t cx e 2 ' ^ eVx v) , which is just Eq. 17 in 
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FIG. 2. Size of the first-order NA transition (t) as a func- 
tion of 10CB concentration (x) in the 8CB-10CB system. 
Comparison between the RG result (full line), HLM (dashed 
line), and experiment. Filled circles are data from Yethiraj 
and Bechhoefer, while hollow circles are latent-heat data from 
Marynissen et al, converted into equivalent size t. 

By tuning the concentration x of 10CB liquid crystals 
in 8CB, one can vary the temperature range of the ne- 
matic phase, and, in effect, tune the parameter b in Eq. 
pi. For larger concentrations, b is negative and the transi- 



tion is more strongly first order. The latent-heat data for 
x > 0.42 are well-fit by mean-field theory (Eq. ||) with 
the parameters c = 1, a = 3.35, b = bo(x* — x), with 
b = 0.395 and x* = 0.42, and e 2 = 0.0421 §|-§4j. For 
x < x* , the quartic term in Eq. (^) becomes positive, 
but the transition remains first order, in agreement with 
the HLM theory. At smaller concentrations, t continues 
to decrease, as expected, but there is a clear deviation 
from the HLM result pj[ . To attempt to fit the data 
in the whole concentration range < x < 0.65, we take 
the above parameters to set the initial values of the cou- 
plings, and then evolve t, e, and b according to Eqs. |(|||, 
with u = v = 1/4, x = 1/16, y = (2V2 + l)/8, z = 1/2, 
and w — l/(2\/2) jlijl - We neglect the change of c un- 
der scaling, since the results are quite insensitive to its 
precise value. The only free parameter left is the final 
value of the temperature to a t which it becomes safe to 
neglect the OP fluctuations. We chose to = 0.01, an or- 
der of magnitude larger than the largest measured t in 
the experiment. The quality of the fit turns out not to 
be critically dependent on this choice. The fit to our the- 
ory and the comparison with HLM is shown on Fig. |^. 
Note that for x > 0.35, there is very little change from 
the HLM result p5| , but as b increases at smaller con- 
centrations, deviations from HLM become significant, in 
agreement with the experiment p6] | . One expects that if 
one could increase b further, eventually t would go to zero 
in accord with Eq. ( pi] ) . The full non-perturbative struc- 
ture of the flow diagram in Fig. [I], with the tricritical and 
the inverted-XY fixed points, however, is beyond the sim- 
ple one-loop (3- functions we used [^Tj. The data suggest 
that, at small charge, the tricritical point is likely to lie at 
some k c > which corresponds approximately to the 
smallest concentration used in the experiment (x = 0). 

Although anisotropy is known to be important in liq- 
uid crystals [^) , its inclusion, besides introducing two new 
couplings, does not essentially change the structure of 
Eqs. 0- H In particular, the sign of the e 4 term in 
Eq. g stays negative, and our main point that the size of 
the transition in good type-I materials is larger than the 
mean-field prediction remains valid. Some small quan- 
titative differences from our result would be expected, 
however. 

In conclusion, we have shown that the naive expec- 
tation that order-parameter fluctuations would always 
decrease the strength of the first-order transition pre- 
dicted by Halpcrin, Lubensky, and Ma is incorrect for 
good type-I materials. While the transition strength docs 
decrease monotonically as the material becomes more 
type II, it becomes smaller than the HLM prediction only 
near the tricritical point that separates type-I and type- 
II regimes. By using one-loop renormalization group to 
calculate the evolution of Landau coefficients in the free 
energy, we can account for the deviations from HLM pre- 
dictions that were observed in a recent experiment. 



Our arguments also lead one to expect a non-trivial 
exponent near the tricritical point where the transition 
eventually should become second order. One would ex- 
pect this point to be accessible experimentally. For ex- 
ample, one can view the variation of concentration x in 
the 8CB-10CB experiments as a way of continuously ad- 
justing the effective molecular length (between 8 and 10 
CH2 lengths). Mixing 8CB with a shorter molecule (for 
example, 6CB ^9|) might then give access to the tricrit- 
ical point and allow at least a crude measurement of the 
crossover exponent <f>. 
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